On 
On 



Route to Chaos in the 15 NH<s far infrared ring laser 

R. Dykstra, M.Y. Li and N.R. Heckenberg 
Department of Physics, 
The University of Queensland, 
St Lucia, QLD. 4072. 



Abstract 



^ ■ The route to chaos as the pump power is raised in the Lorenz-like 15 NH3 laser 

■ is studied using transients and compared with the route which is expected in the 

. Lorenz Haken equations. The differences between the routes to chaos, as normally 

displayed by the experimental system and the Lorenz Haken equations, is described 
I in terms of the bifurcations experienced in each system. It is also shown that the 

experimental system has a small parameter regime where the route to chaos is the 

same as in the Lorenz Haken equations. 
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O ! 1. INTRODUCTION 

In 1963 LorenzS reduced the equations for convective fluid flow into three first order coupled 
nonlinear differential equations and demonstrated with these the idea of sensitive dependence 
upon initial conditions and chaos. It was shown by Haken in 1975 that these equations are 
isomorphic with the Maxwell Bloch laser equations for a homogeneously broadened single 
mode travelling wave resonantly tuned laser in which the modes are presumed plane waves 
and the pumping is presumed uniform!. Even though these equations are isomorphic their 
conventional layout is different, as shown in equations [T]-|3] and |J-^. 

x = a[y-x] (1) 
y = px — y — xz (2) 
z = xy — (3z (3) 

E = kP-kE (4) 
P = 1± [ED-P] (5) 
D = 1{l [\-1-D-XEP] (6) 

The variables in the Lorenz equations, namely x,y and z correspond to the slowly varying 
amplitudes of the electric field E and polarization P and the inversion D respectively in the 
Lorenz-Haken equations. The parameters are related via /3 = a = and p = A + 1, 
where 711 is the relaxation rate of the inversion, 7^ is the relaxation rate of the polarization, k 
is the field relaxation rate, and A represents the pump power, normalized to the threshold 
for CW lasing. 
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It was conjectured that Lorenz-like chaotic behaviour could exist in an optically pumped 
far infrared ring laser. Subsequently, such chaos, characterised by variable length trains 
of pulses of increasing amplitude ("spirals") was observed by Weiss and coworkers in NH 3 
lasers, under appropriate conditions!! In this paper we will not discuss the other sorts of 
chaos exhibited by the laser while it is under other conditions. Detailed study has since shown 
that there are only a few minor differences between the laser's Lorenz-like chaos and Lorenz 
Haken chaos@. One of the most obvious differences is the appearance of an anomalously 
large peak nearly always observed at the start of every spiral in the laseill. This has been 
attributed to deviations from the simple model due to Doppler broadening or coherent 
pumping effects on the pump transition!. It was also found that this effect appears if the 
theory simply accounts for the fact that the laser and the pump have transverse structures 
(i.e. there exist real modes in the cavity) and are not plane wave as the Lorenz Haken 
equations assumed. Perhaps it is a generic effect of any small deviation from the Lorenz 
equations. The other obvious difference in the behaviour of the laser is that under nearly all 
circumstances the laser pulses periodically below chaos threshold and does not make a direct 
transition from the steady state to chaos as the Lorenz Haken equations predict. Although 
topological arguments suggest that in any continuous transition from a fixed point to a 
strange attractor the attractor should become one-dimensional, corresponding to a periodic 
orbit, at some point, the transition to chaos in the Lorenz equations is a global bifurcation 
with no intermediate states. Further, the model accounting for nonuniformity of pump 
and electric field also shows a transition with no intermediate periodic state. The aim of 
this paper is to describe the route to Lorenz like chaos in an optically pumped far infrared 
ammonia laser and suggest the possible nature of the bifurcations and hence the underlying 
topology of the laser attractor as the pump power is increased from the lasing threshold to 
chaos. 

The next section will give a short revision of the route to chaos with increasing pump 
power for the Lorenz Haken equations to allow easy comparison with the route to chaos 
observed in the laser, which will be described in section |3]. In section ^ the results of section 
^| will be extended to include the effect of detuning in the laser. 



2. THE LORENZ HAKEN EQUATIONS 

The Lorenz equations have been studied extensively over the years as a basic illustrator 
of chaos and demonstration platform for various theories relating to chaos. An extensive 
review of Lorenz chaos was undertaken by Sparrow!, who subsequently showed numerically 
that there is considerable complexity in the dynamics of the Lorenz equations which was not 
previously analytically explained. This included such phenomena as homoclinic explosions 
which give rise to stable periodic motion, and a heteroclinic bifurcation. 

To explain how the route to chaos in the laser differs from that of the equations a 
brief revision of the bifurcations in the equations will be given. Below lasing threshold the 
origin is globally attracting. That is, all trajectories finish at the origin ifO<A+l<l. At 

A > a simple bifurcation occurs which creates two stationary points at (±y^j^, ±y^^, A) 

henceforth called C\ and Ci- At this bifurcation the origin loses its stability, however C\ and 
C*2 are stable. The two points represent the steady state solutions for the system in the two 



2 



possible orientations that they might find themselves, with the elctric field and polarization 
differing by a /Pi phase shift. At A + 1 = -^-(-^z^z^ - ) the equations go through a second 
bifurcation in which the points C\ and Ci lose their stability. The bifurcation is of the type 
known as a Hopf bifurcation, which may occur in two ways. The bifurcation is "supercritical" 
if each point loses its stability by expelling a stable periodic orbit. It is "subcritical" if they 
lose their stability by absorbing an unstable periodic orbit. For values of A greater than the 
Hopf bifurcation point the origin and both C% and C<i are all unstable. There is another 
point on the route to chaos which delineates an important event and this is the occurence 
of a homoclinic orbit. The homoclinic orbit separates the region between normal stable 
behaviour and possible metastable chaos, associated with an unstable periodic orbit. This 
can be understood by considering operation just above lasing threshold. A trajectory started 
very near to the origin will move out and, as shown in figure 1, quickly spiral into either 
C\ or C2 depending on which the point started closer to. As A is increased the trajectories 
around the Cj widen until the trajectory in both the forward and backward propagation of 
time moves towards the origin. There is therefore a homoclinic orbit associated with the 
origin. If A is increased slightly further the trajectory jumps across to spiral into the opposite 
stable point. Figure 1 illustrates this point and SparrowH gives a more detailed explanation 
of how this occurs. Increasing A above the homoclinic explosion, as shown by Sparrow, 
the homoclinic orbit becomes an unstable periodic orbit, getting smaller and finally getting 
absorbed by the stable point C\ or C2 at the Hopf bifurcation. 

The implications of the existence of these orbits is best described by figure 2 where the 
bifurcations are placed into context with each other. When the unstable periodic orbits 
decrease in size there is a possibilty for the trajectory to oscillate between the two sides 
of the attractor and do so for long periods of time. However because there are two stable 
attractors embedded in an unstable attractor the trajectory must finally intersect the stable 
manifolds and spiral into the stable points C\ or Ci- So instead of the infinitely long single 
trajectory as experienced in chaos there are an infinite number of finite length trajectories. 
This is generally called "preturbulence" or "metastable chaos" and refers to the fact that 
the dynamics may look chaotic, but will at some finite time intersect the stable manifold 
and spiral toward one of the stable points. 



3. ROUTE TO CHAOS IN THE LASER 

The experimental setup we use has been fully described by Tin Win et afl Briefly, it is 
a 15 Ai?3 far infrared ring laser which is pumped by a 13 C02 laser whose pump output is 
attenuated by an acousto optic modulator (AOM) of rise time 0.8 /j,s. To gain more insight 
into the attractor properties, the laser was switched from above chaos threshold to below 
chaos threshold periodically by the AOM. Each period above chaos was about 50-100 chaotic 
pulses long, whereas the period below chaos threshold was generally about 100-150 pulses 
in length. The time between pulses is about 0.8 fis. Figure 3 illustrates the time series 
of the Lorenz Haken equations as a result of such a sudden decrease in A from just above 
chaos threshold to below chaos threshold. The transition to the new state involves two 
phases. Firstly there is a variable length period of chaotic pulsations, followed by a damped 
oscillation. Theoretical studies of this behaviour0 have shown that the length of the chaotic 
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phase can be characterised by an exponential distribution with a mean length dependent 



shows a typical time series for the 15 NH% far infrared laser, demonstrating the predicted 
behaviour. It should be noted that the only variable which is measurable in the laser is the 
intensity. This of course places some restrictions on the information which is processable 



shown that the chaos exhibited by the laser, for specific parameters is very similar to Lorenz 
chaos. Because the chaos is Lorenz-like it can be assumed that the dynamics above chaos 
threshold must have a topologically similar structure to the Lorenz Haken equations. At the 
other extreme, Tin Win et aM determined the parameters 7y, 7_i_ and k for the 15 NH3 laser 
by comparing the decay envelopes for transients for low pump power with the linearised 
Lorenz Haken equations. These studies showed that the laser operates in reasonably close 
agreement with the Lorenz Haken equations for low pump powers. 

Figure 3 and figure 4 compare the dynamics of the experiment with what is expected from 
the equations. In the parameter range where the dynamics of the laser exhibits preturbulence 
its behaviour also seems to be very similar to the Lorenz Haken equations^. In the Lorenz 
Haken equations, close to but above the homoclinic explosion a trajectory started in the 
vicinity of the origin should go around one of the points Cj once and then spiral into the 
other point. It is normally possible to measure only the intensity in these experiments which 
obviously makes it difficult to observe a jump from one side of the attractor to the other. 
One way to detect a jump is to measure the phase of the laser radiation by heterodyne 
detection. We have carried out this experiment using a second 15 NHz laser (a standing 
wave laser operated continuously with a few milliwatts of output power) as a local oscillator 
and a Schottky barrier detector as mixer with the two beams combined at a polyethylene 
beamsplitter. The output from the detector contains a homodyne spectrum from the ring 
laser, extending out several megahertz, and the heterodyne spectrum centred on the mean 
frequency difference between the two lasers which was limited to about 6 MHz. The detector 
signal was digitised, the homodyne spectrum removed numerically and quadrature field 
components (or amplitude and phase) reconstructed by numerically mixing with a signal at 
the mean difference frequency, in analogy with RF engineering methods. Appendix A has a 
thorough explanation of the method used to extract the amplitude and phase of the signal. 
This method is similar to that used by Tang et al0, but handles transients better. Examples 
of the results are shown in figures 5 and 6. In both cases the laser pump power was rapidly 
switched on from zero to a level giving CW output, and the transient approach to the steady 
state was recorded. Figure 5 shows the intensity and phase variation at low pump powers 
with an apparent slow drift in phase (a result of a drift in local oscillator frequency relative 
to the ring laser) but no jump in phase. Figure 6, however shows the intensity and phase 
variation with the pump power above the first homoclinic explosion, with a 7r phase jump 
clearly present between the first and second pulses. 

The only difference between the laser and the equations is the occurence of a stable 
oscillation within the preturbulent regime in the laser. Instead of finally decaying to a steady 
value corresponding to a fixed point, the decay is to a stable oscillation. This oscillation 
is also met if the pump power is raised slowly in a quasi-static manner, with a threshold 
somewhat below the chaos threshold. The oscillation starts small at lower pump powers 




from any time series 



obtained. Comparison of return maps and statistical analysis^ have 
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and grows with increasing pump power until the laser shows no more preturbulence. When 
the stable oscillation grows to be about the same size as the preturbulence, the laser jumps 
into chaos and stays there. Figure 7 shows a typical graph of percentage pump power below 
chaos threshold and oscillation amplitude for a particular set of parameters. Figure 4 shows 
an example of such behaviour. When this is compared with the time series of figure 3 it is 
obvious that this does not occur in the Lorenz Haken equations. 
To revise what happens as the pump power is increased; 

1. For pump powers below 120 mW at a pressure of 38 \ibar there is no lasing 

2. For powers greater than this the laser makes a transient oscillation which damps 
to steady state much like the Lorenz Haken equations, moving towards one or other of the 
stable fixed points (Figure 5). 

3. A point started near the origin and close to Lorenz- like chaos will exhibit a jump 
from one side of the attractor to the other and subsequently collapse into steady state lasing 
(Figure 6). 

4. For slightly higher pump powers the laser exhibits preturbulence much like that of 
the Lorenz Haken equations. The differences between them lies in the appearance of a stable 
oscillation in the laser which grows in amplitude as the pump power increases (Figures 3 
and 4). 

So far the dynamics of the Lorenz Haken equations and the dynamics of the laser have 
been described to explain their respective routes to chaos as the pump power is increased. 
Both routes look quite similar except for the appearance of the stable oscillation in the laser. 
A possible explanation for the occurrence of such a stable oscillation is shown in figure 8. 
Starting at low pump power the first significant event, as pump power is increased, is the 
appearance of the two unstable periodic orbits around the stable points Cj as is the case with 
the Lorenz Haken equations. Increasing the pump power further the points Cj bifurcate in a 
supercritical Hopf bifurcation, each expelling a stable periodic orbit and becoming unstable. 
A representative diagram showing all the orbits for a specific pump power is shown in figure 9. 
Increasing the pump power from here on results in an increasing size of the stable homoclinic 
orbits. Chaos is reached when the unstable and stable homoclinic orbits annihilate each other 
to leave Lorenz like chaos. The most trivial implication of this route to chaos is that the 
points Ci go unstable at lower pump powers than they do in the Lorenz Haken equations. 
This also has the implication that the range of pump powers over which preturbulence exists 
may be reduced. 

4. THE EFFECTS OF DETUNING 

During the course of the experiment it was observed that the amplitude of the stable os- 
cillations during preturbulent chaos is dependent upon detuning of the laser cavity. It was 
found that within an extremely small range of detuning, preturbulence occurred with the 
laser relaxing down to steady state rather than the stable oscillations. This range was far 
narrower than the range over which Lorenz like chaos was observed, and so narrow that 
it was not possible to carry out an exhaustive investigation. However, figure 10 shows an 
intensity time series of such an event. At the moment it is not possible with our setup to 
measure how far from zero detuning the laser is, no conclusions can therefore be made as 
to where this region resides. It would not be surprising however to find this region at zero 
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detuning. This would therefore mean that there does exist a regime where the laser chaos 
corresponds exactly to the Lorenz laser case described by the Lorenz Haken equations. Such 
a regime has also been observed in the U NH 3 laserEl. 



5. CONCLUSIONS 

It has been shown, based on measurements of the transient behaviour of the laser, that the 
15 NH 3 far infrared laser has a route to chaos with increasing pump power which is slightly 
different from that of the Lorenz Haken equations. As with any other laser it has a minimum 
pump power below which no lasing is possible. For powers greater than this the laser makes 
a transient oscillation which damps to steady state much like the Lorenz Haken equations. 
Using heterodyne techniques to see the field there was found to be a pump power above 
which the laser exhibits a jump from one side of the attractor to the other where the field 
subsequently collapsed to steady state lasing. At slightly higher pump powers the laser was 
found to exhibit preturbulence leading to a stable oscillation which grows in amplitude with 
increasing pump power. Finally, above a certain power Lorenz like chaos is observed. 

The route to chaos, as observed in our laser, has been conjectured to involve an extra 
stable periodic orbit which annihilates with the unstable periodic orbit to leave Lorenz like 
chaos. 

It was also shown that for an extremely small range of detunings the laser will show 
preturbulence without the occurence of the stable oscillations. These phenomena will provide 
an excellent test for improved theoretical models of the laser. 



6. Appendix A 

In the heterodyne measurement of the chaotic laser, the signal from the chaotic laser and the 
signal from the local oscillator are coupled together by a beam splitter. The mixed signals 
are then detected by a Schottky barrier diode. The electric field of the chaotic laser and the 
local oscillator are: 

E x {t) = A x {t) wo[u x t + (pit)} (7) 
E 2 (t) = A 2 sin(u 2 t) (8) 

where A x , uj x and are the amplitude, the "optical" frequency and the phase of the electrical 
field of the ring laser respectively. A 2 and u 2 are the amplitude and the "optical" frequency 
of the local oscillator respectively. A x and 6 are time dependent. A 2 , uj x and u 2 are constant. 
The response of the Schottky diode isEj 

I{t) oc e aV = e ^ +E ^ = 1 + P(E X + E 2 ) + ^ 2 {E X + E 2 f + ... (9) 
where a and (3 are constants. 

Because the rest of the electronics responds only to frequencies much lower than the 
"optical" frequencies u x and u 2 , and a bias voltage is applied to the Schottky diode to make 
it favor the square term, in the output from the Schottky detector only the square term 
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is important. So for At ^> l/u>i, l/a>2 (At is the response time of the Schottky diode) the 
output from the Schottky diode is 

< I > At cc < (E 1 + E 2 f > 

= < A\ sin 2 (cuxt + <p) + A\ sin 2 (^ 2 t) + 2AiA 2 sin^t + <p) sin(u 2 t) > 



The output from the Schottky diode is a mixture of the desired heterodyne signal carried 
by the beat between the chaotic laser and the local oscillator with an unwanted homodyne 
component. The time variation of A% and lead to a finite width for both the homodyne and 
heterodyne spectra in the frequency domain. It was ensured that the spectra did not overlap 
so that the homodyne and heterodyne signals could be separated by using simple filters. If 
they overlapped a more complicated detection setup, for instance using two Schottky diodes 
to detect the homodyne signal and the mixed signal respectively, would be necessary. 

In the heterodyne signal, A\ and <fi are carried by the beat between the chaotic laser and 
the local oscillator. The carrier signal must be removed from the heterodyne signal in order 
to reconstruct the chaotic laser field. A Fast Fourier Transformation (FFT) was applied to 
the experimental record to obtain the signal spectrum and determine the cutoff frequencies 
for a low and high frequency filter. These filters were subsequently used on the original 
experimental record to separate the homodyne and heterodyne signals. The homodyne signal 
was constructed directly from the low pass filtered portion of the signal. The heterodyne 
time series was separately mixed with a sinusoidal and a cosinusoidal signal of the beat 
frequency and low pass filtered in order to reconstruct the two orthogonal components of 
the heterodyne time series and hence the field components and phase. 

Tests on synthetic data show this method is reliable for in analyzing transient signals, 
as it is able to reconstruct the original field amplitude and phase perfectly. By contrast, the 
signal processing method used by Tang et al§, which is based on manipulation of Fourier 
transforms, is able to reconstruct periodic data reasonably well, but it introduces errors in 
the amplitude and the phase when the signals are transient. 




(10) 
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Figure Captions 



Figure 1: Trajectories of the Lorenz Haken laser equations at pump powers, A equal to 1.0 
and 2.0 below the first homoclinic explosion, ^3.0351841193 for the homoclinic orbit, and 
4.0 above the homoclinic explosion. Parameter values of 711 = 1/8 , 7± = 1/2 and k — 1 are 
appropriate for the 15 NH 3 laser and give a chaos threshold of 13. The circles indicate the 
stable points towards which the trajectories will gravitate. 

Figure 2: The bifurcation diagram for the Lorenz Haken laser equations. The x = y axis 
symbolises the symmetry along that plane of the equations. 

Figure 3: Intensity time trace of preturbulence in the Lorenz Haken laser equations at 
7|| = 1/8 , 7j_ = 1/2 and k — 1 resulting from the fast switch down from A = 14 above the 
chaos threshold to A = 10 below the chaos threshold. 

Figure 4: Typical intensity time trace of preturbulence in the 15 NH 3 laser as the laser is 
switched from slightly above the 5.85 W chaos threshold to a factor of 0.1 below chaos 
threshold. 

Figure 5: Phase and intensity evolution of the 15 NH 3 laser. The dashed line shows the 
decaying pulses of the intensity and the solid line the phase of the light. The phase shows a 
slow drift as the pulses decay indicating the absence of a jump from one side of the attractor 
ot the other. The pump power was a factor of 0.7 below chaos threshold. 

Figure 6: Phase and intensity evolution of the 15 NH 3 laser at switch on. The dashed line 
shows the decaying pulses of the intensity and the solid line the phase of the light. The 
phase shows a distinctive jump of it after the first pulse to indicate a jump from one side of 
the attractor to the other. The pump power was a factor of 0.5 below chaos threshold. 

Figure 7: Amplitude of the intensity oscillations observed in the 15 NH 3 laser as a function 
of pump power. 

Figure 8: The bifurcation diagram of the l5 NH 3 laser. The difference between this diagram 
and that for the Lorenz Haken equations shown in figure 2 is the presence of a supercritical 
Hopf bifurcation within, and which finally annihilates, the already existing homoclinic orbit. 

Figure 9: Representation of the stable and unstable orbits relative to their respective fixed 
points. The format of the figure is as for figure 2. 

Figure 10: Preturbulence in the the l5 NH 3 laser for the detuning range where the pretur- 
bulence settles down to a steady state. The laser is switched from slightly above the 5.85 
W chaos threshold to a factor of 0.05 below chaos threshold. 
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